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Abstract 

Boundary conditions for a superconducting order parameter at a diffusive 
scattering boundary are derived from microscopic theory. The results indicate 
that for all but isotropic gap functions the diffusive boundary almost com- 
pletely suppresses surface superconductivity in the Ginzburg-Landau regime. 
This indicates that in anisotropic superconductors surface superconductiv- 
ity can only be observed for surface normals along high symmetry directions 

where atomically clean surfaces can be cleaved. 
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Superconductivity in UPt 3 is believed to be described by an unconventional order pa- 
rameter which has two complex components rji and rj 2 [1-5]. This description allows for an 
explanation of the many unique experimental features associated with superconductivity in 
this hexagonal material. For example, the upper critical field for fields in the basal plane 
(H^) as a function of temperature displays a kink at a temperature T* [6]. At temperatures 
below T* one component of the order parameter orders at H^, while for temperatures above 
T* the other component orders at H^K The observation and experimental investigation [6-8] 
of surface superconductivity in UPt 3 has lead to an examination of this phenomenon in un- 
conventional superconductors [9-11]. In contrast to , the upper critical field for surface 
superconductivity for fields in the basal plane (H^) exhibits no anomaly with tempera- 
ture [8]. This can be understood if one component of the order parameter is suppressed 
at the surface, allowing surface superconductivity to occur only with the other component 
[9,11]. It is therefore important to understand under what conditions the superconducting 
order parameter is suppressed at a boundary. Microscopic calculations for boundary con- 
ditions at specular reflecting surfaces have been conducted [10,12]. It is also interesting to 
examine the effects of diffusive boundaries on surface superconductivity. It is known that 
rough surfaces are pair breaking for anisotropic superconductors [13,14] and therefore are 
expected to suppress surface superconductivity. Here we investigate how strong such a sup- 
pression will be for general order parameter symmetries. First, we examine the solution of 
the isotropic Ginzburg-Landau model with general boundary conditions to gain an under- 
standing of the effects of the boundary conditions on surface superconductivity. Then we 
use a weak coupling microscopic theory to calculate the boundary conditions at a diffusive 
scattering boundary for general order parameter symmetries. 

The Ginzburg-Landau free energy density for an isotropic superconductor in an external 
magnetic field H is 

F = a (T - T c )|Vf + K(Dip) ■ (Dip)* + h 2 /8vr - h • H/4vr (1) 

where D = d — (z2e//ic)A, and h = d X A. We have only kept terms to order if) 2 since 
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we are interested in determining the upper critical field. In the presence of a surface the 
following surface free energy density is added 

Fs = gW- (2) 

The Ginzburg-Landau equations are 

atp = kB ■ Dtp h = H (3) 

where a = a (T — T c ) with the boundary condition n • Dif)\ sur j ace = (l/b)ip\ sur f ace where 
n is the surface normal and b = n/g is the extrapolation length and is described in Fig. 1 
(also see Ref. [18]). For an applied magnetic field orthogonal the the surface normal this 
equation can be solved by following the method of Saint James [15]. The solution is 

^ = l - (4) 

where H C2 = 27r</> /<f; 2 (T), <f> is the elementary flux quantum, £(T) = \J —n/ a is the coherence 
length, and //(/) is given by 

/°°(2u -/i- /)e-("-^V( 1+ ^ 2 -' 2 )/ 2 du = 0. (5) 
Jo 

where / = (H C2 / H C3 )^ 2 C(T) /b. The ratio H C JH C2 function of temperature is given in 
Fig. 2. 

Microscopic calculations are required to determine the ratio b/£(0). Such calculations 
have been performed for general order parameter symmetries in the presence of a specular 
reflecting surface in a weak coupling model with a spherical Fermi surface [10,12]. Here 
similar calculations are performed for a diffusive scattering surface (this has previously been 
done for isotropic [16], weakly anisotropic [17], and p-w&ve [13] order parameters). The 
correlation function method developed by deGennes [18,13] and extended to unconventional 
superconductors by Sigrist and Ueda [19] is used. We consider a weak coupling model with 
a spherical Fermi surface and assume that there is no spin flip scattering at the surface. The 
development of the formalism initially parallels that of Sigrist and Ueda [19]. 
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The Hamiltonian i 



is 



H=J2(ks\H \k's')cl s C k , s ,+ J2 ^1^2, S 3, S 4( k , k/ ) C i/2-k, Sl C q/2 + k, S2 C q/2 + k', S3 C q/2-k', S4 ( 6 ) 



k ,k ! k,k ! ,q 



s l. s 2. s 3. s 4 



where H is the single particle Hamiltonian including the interaction due to the 
boundary and ct jS destroys a free electron with Momentum k and spin s. Defining 
<?,,,,(k,k',T) = -(r T {c k , s (r),4, s ,(0)}) and i^,(k, k', r) = (T T {c{,, ,„ (r)4,,(0)}) where 
c k,s(''") = exp Hr Ck, s exp~ Hr and T T is the imaginary time ordering operator and finding the 
equations of motion for these averages leads to the linearized gap equation 



A S1 , S2 (q, k) =-k B Tj2 E Vs 2 , Sl ,s 3 , Si (K k')Gl S5 (k' + q/2, k" + q'/2; 



Un k',k",q' 

s 3. s 4. s 5. s 6 

^0 



xG° 4jS6 (-k' + q/2, -k" + q'/2; -*a; n ) A, 8i , 6 (k", q') (7) 

where u n = TrkBT(2n + 1) are the Matsubara frequencies and the normal state electron 
Greens function G° s s ;(k, k', iu n ) is given by the Fourier transform of 

G SjS ,(R,R,«^) = 2^ : S s,s> (8) 

where </>j, jS (R) are the eigenfunctions of the single particle Hamiltonian Ho- We assume that 
the interaction can be written in the weak coupling form 

y sl , S2 , S3 , S4 (k,k') = E 5 (r)A S2iSi (r,m,k F )At 3iS4 (r,m,k' F ) (9) 

r,m 

where e(k) and e(k') are restricted to lie within an energy e c of ep, T refers to an ir- 
reducible representation of the point group and m to the basis of the Y representa- 
tion. If only one representation is important then the gap matrix can be written as 
A(R, kp) = J2 m ?7m(R-) A m (kp). Fourier transforming the linearized gap equation with 
respect to the center of mass variables q and q' , substituting the above forms for the 
gap function and the potential, and using the following orthogonality condition for the gap 
matrix 

—5-2- / J 2 A ; rr[A / (k r )A m (k r )] = 28 m<l (10) 



where the integral is over the Fermi surface, gives the following equation for the order 
parameter r)i(R) 

r)i(R) = J d 3 R'K tl (TL } TL') Vl (TL') (11) 
with A% (R,R') given by [19] 

A%(R,R')= dt dei^ (12) 

? (H <rA « T [^ J(R) ] Aj S J(R/ ' t} l I 17 ) 15(6 " e ^ (13) 

where J(R) is the current operator. We use the semiclassical and weak coupling approxi- 
mations, which entail Y^{trk\k 3 )6{e v - e) « N{{)){trk\k 3 ) t=tF , classical [18,13] where N(0) 
is the density of states at the Fermi surface, to arrive the following from for the kernel [19] 

K l3 =gN(^k B TY: r^exp- 2 l-l*(trAt[^J(R)]A4^J(R',t)]\ (14) 

ujn J ° \ 2k F 1 L 2k F 1 I e F , classical 

where the expectation value is an average in a canonical ensemble for an electron with 
momentum on the Fermi surface. In the presence of a single boundary, the kernel has two 
parts, a direct contribution (K d ) which is the contribution when no boundary is present and 
a contribution due solely to the scattering at the boundary (K r )- For simplicity, we assume 
a spherical Fermi surface in which case the direct contribution is given by Sigrist and Ueda 
to be [19] 

A £( R ) = — ^ — E — — jr 2 - ex p VF ■ ( 15 ) 

The transition temperature is given by the condition / ePi^fQ-(R) = 1 where i corresponds 
to only one component of the order parameter (for a more detailed discussion of this point 
see [18]). For a diffusive boundary the expectation value in K r - is given by 

73^ , m l 



trA] A,) = (4iV(0)4)- 1 / 7^4 / ^V^ 3 [ R/ - R± + P 1 --P'( 1 + -)] 

' J ( Ztt ) J wp p z m p z 



Pz<0 p' z >0 

6(p - p')S(f- - e F )tr [AJ(p) A,(p')] • (16) 
Am 1 J 
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In this equation a quasiparticle that had initial position R = (R± } R Z ) (z is the component 
along the surface normal) and initial momentum p = (p±,p z ) has been scattered by the sur- 
face emerging with momentum p'. The p' z l^p' represents the probability of emerging from 
the surface with momentum p' (note this is independent of p) and the 3D delta function 
gives the time dependent position R'(t) of the quasiparticle given that at time t = —mR z /p z 
the quasiparticle is at the surface with position Rj_ — p±R z /p z (these correspond to the po- 
sition on the surface and the time required to reach the surface for a quasiparticle with 
initial position and momentum given by R and p)(see Ref. [16] for a discussion of diffusive 
scattering for the case of isotropic superconductors). We assume that the order parame- 
ter varies only in the direction along the surface normal and therefore wish to determine 
/ d 2 (R ± - R' 1 )K lJ (R^') = Kij(z,z'). The resulting kernel is 

v f s A\(s)\ x / dsexp v f s 



+7r 1 t r [ j dsexp v f s A-(s)J x y j dsexp v f s Aj(— s)J j (17) 

where F tJ (s) = (l/2)tr J 2v d<j> [AJ(s, <f>)kj(s, <j>) + k\{-s,<j>) k 3 (s^% k t (s) = 
Jq V d(f>Ai(s } </>), and A;(</>, s) is given by setting k = (\A — s 2 cos <f>, \J\ — s 2 sin </>, s) in A(k). 
To obtain Eq. 17 the following was used 

[ 2V d<f>k](<f>, p, z)k 3 (<i>, p, z) = r d<f>k](<f>, p, -z)k 3 (<i>, p, -z) (18) 
Jo Jo 

where A((f> } p } z) is given by setting k = (p cos </>, p sin </>, z) in A(k). Eq. 18 arises because 
both A; and A 3 transform identically under parity. A similar development for a specularly 
reflecting surface gives the same result as Samokhin [10]. 

To proceed we consider the i = j contributions only (these are frequently the only 
contributions along high symmetry directions). After introducing x = z/^ 0} where £ = 
HvF/2TrkbT C} the integral equation for the order parameter becomes 

poo „ 

rji(x) = / dx'K(x, x')rn(x') (19) 



where K = ifivp /2tv kbT c )Ka. It can be verified that rji = rji (l + x^ /bi) is a solution to Eq. 
19 as x — > oo. As pointed out by deGennes [18], the linear dependence of the order parameter 
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on x appears to give an unphysical result as x — > oo, however nonlinear terms neglected 
in Eq. 19 will introduce a negative curvature to the order parameter so that it will achieve 
its bulk value for ^ x > £(T). We are interested in the region ^ x ~ Co << C(^) where 
this curvature is negligible. To find the coefficient 6;/Co we use the variational approach of 
Svidzinskii [20] as it is presented by Samokhin [10] and by Barash et.al. [12]. Substituting 
i) = C(x + q(x)) (then 6/Co = lim^oo q(x)) into the Eq. 19 gives 



E( 



x 



oo 



q(x) = h / K(x,x')q(x')dx' (20) 

2 Jo 

with E(x) = 2 J °° x'K(x,x')dx' — 2x. The above equation can be found by minimizing the 
functional 

'!'[ 1 = ^ dx( ii x )\- ( ii x ) ~ Io° dx'K(x, x')q(x')] 
Wl ~ [/ °° dxq{x)E{x)Y ' 1 ' 

The minimum value of ^[q] is given by 

* mm = 2j oo dxq(x)E(xy (22) 



The coefficient b can be expressed in terms of ^ m in as 



1 Jq 00 dxxE(x) + 



2 JO i 4 ^ 



Co \ Jq 00 dxE(x) — J'q 00 dx'x'[j^ dxK(x', x) — 1] 



(23) 



Using a constant for q(x) gives the following result 
bi (7C(3))" 1 



Co ft s*F(s)ds + ±tr [/o 1 sA](s) x £ s^A t (-s) - J, 1 s*A}(s) x £ sA t (-s) 

x I — \ T s 3 F(s)ds + n-Hr f 1 s 2 A](s) x C s 2 A t (-i 
24 Uo Jo Jo 



(24) 



(7((3)) 2 [lo s 2 F(s)ds + ir-Hr J, 1 sA](s) x £ s 2 A t (-s) 
+ 2** |i sF(s)ds - Tr-Hr J, 1 sA}(s) x £ sA t (-s)] 

where ((3) = £ n l/(2n + l) 3 . 

Values for 6;/Co are given in Table 1 for various functions A(k) corresponding to irre- 
ducible representations of Deh for a surface normal along the hexagonal a direction. For the 
case of p-wave pairing an exact solution can be compared to the variational solution. The 
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p-wave order parameter transforms as a vector under spacial rotations. For the order param- 
eter component transverse to the surface normal the variational solution gives b t /^ = 0.53 
which compares favorably to the exact result b t /^ = 0.54 [13]. For the longitudinal compo- 
nent the order parameter obeys r)i(0) = irrespective of the form of the boundary [13]. In 
this case the variational approach gives &;/£o = 0.11 which is non-zero. This non-zero value 
arises because rj = C(b/^ + x) is an asymptotic solution to Eq. 19 and the exact boundary 
condition is valid only on the surface (see Fig. 1). The variational result for the constant 
order parameter (&/£o = oo) is exact. 

Note that in general £(0) = \J 'k/ 'a T c ^ £ (for isotropic superconductors £(0) ~ 0.2^ o 
[18]). However £(0) ~ £o and the values of b/£(0) (~ 0.6) give rise to H c? H C2 to less 
than a tenth of a percent within the temperature range shown in Fig. 2. This indicates 
that diffusive scattering effectively completely suppresses surface superconductivity for all 
but isotropic order parameters. Since the electronic wavelength is typically on the order 
of atomic length scales a surface will usually be diffusive. An interesting implication is 
that surface superconductivity is expected only to occur on surfaces with normals along 
high symmetry directions, where atomically clean surfaces can be cleaved. These results 
are consistent with the observation by Keller et. al. that cutting the crystals destroyed 
surface superconductivity in UPt 3 [8]. It has been proposed that turning a cylindrical 
superconductor in a magnetic field orthogonal to the axis of symmetry and measuring the 
surface superconductivity can determine the symmetry of the order parameter [10]. Our 
results indicate that such an experiment is not feasible because the surface of the sample 
will be diffusive. 

We acknowledge the support of the Natural Sciences and Engineering Research Council 
of Canada. We also thank N. Keller, J.L. Tholence, A. Huxley, and J. Flouquet for making 
their work available to us prior to publication. 
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TABLES 



r 


A(k) 


— specular 


— diffusive 




iOy 


00 


00 




i(Jy{kl + k 2 y ) 


00 


12 




iuykl 


00 


1.4 


A 2g 


i(Jy(kg — 3k X ky)(ky — 3kgky) 





0.72 


Big 


z(jyk z (ky 3kyk x ^ 


00 


0.66 


B 2 g 


i(Tyk z (k x — 3k x ky) 





0.63 


Elg 


%Oyk z k x 





0.64 




% Q ' y f\t g f\t y 


00 


0.46 


E 2 g 


% £T ' y 2 f\t f\t y 





0.64 




io y (k x — ky) 


00 


1.3 


Aiu 


ia y a ■ zk z 


00 


0.53 




idyCf ■ zk z (k^. — Skgky^ky — 3k 2 ky^j 





0.68 


Bin 


%u yCT ' z(ky t^kyk^^ 


00 


0.72 


B 2u 


10 y(J ' Z^k g t^kgky^ 





0.68 


Eiu 


ia y a ■ zk x 





0.11 




ia y a ■ zk y 


00 


0.53 


E 2u 


%(j y& ' z^k-gkykg 





0.58 




icjycr ■ z(k x — ky)k z 


00 


0.67 


TABLE I. Boundary conditions for gap matrices transforming as selected basis functions of the 
point group D&h. The surface normal is along the hexagonal a direction. The specular boundary 
condition is found by applying P x A(k) = A(k — 2k • xx). If P x A(k) = A(k) then b = oo and if 
P x A(k) = — A(k) then 6 = 0. The 0{ are the Pauli matrices. 
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FIGURES 

FIG. 1. Schematic representation of the extrapolation length b. (a) Microscopic depiction of 
the spacial variation of the superconducting order parameter near a superconductor to insulator 
boundary, (b) Macroscopic representation of (a). 

FIG. 2. H C? J H C2 as a function of reduced temperature for various values of &/£(0). 
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